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Abstract 

We present a new way to tune the electron-phonon couplmg (EPC) in graphene by changing 
the deformation potential with electron/hole doping. We show the EPC for highest optical branch 
at the high symmetry point K, acquires a strong dependency on the doping level due to electron- 
electron correlation not accounted in mean-field approaches. Such a dependency influences the 
dispersion (with respect to the laser energy) of the Raman D and 2D lines and the splitting of the 
2D peak in multi-layer graphene. Finally this doping dependence opens the possibility to construct 
tunable electronic devices through the external control of the EPC. 

PACS numbers: 71.15.Mb, 63.20.Kr, 78.30.Na, 81.05.Uw 
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A large amount of work envisioning exciting new applications of graphene-based devices in 
nanoelectronics has been published in the last years (see Ref. Q and references therein). The 
performance of those electro-optical graphene-based-devices^ is governed to a large extent 
by the electron-phonon coupling (EPC), or, more precisely, by the deformation potential. 
For example in high-current transport the scattering with phonons increases the differential 
resistance in carbon nanotubes and graphene^i"-. It has been widely assumed that the de- 
formation potential in graphene is constant with respect to the electron/hole concentration. 
For a proper description of device performances it is necessary to control the validity of 
this approximation. Indeed, here we show that the deformation potential displays a rather 
strong doping dependence which should be taken into account in the design of new graphene 
devices. 

The interaction between electrons and phonons in graphene and graphite has been 
studied with many experimental techniques ranging from Angle Resolved Photoemission 
Spectroscopy (ARPES)^, inelastic x-ray scattering(IXS)^, Scanning Tunnehng Spectroscopy 
(STS)^ to Raman spectroscopy^. In particular, Raman spectroscopy is commonly employed 
in graphene characterization because it is sensitive to the number of lajersr^'^, the doping 
level^"— and the graphene edges^^. In order to interpret all above mentioned experiments, 
a complete knowledge of the electronic structure, the phonon dispersion and the electron- 
phonon interaction is required. In graphene the electron-phonon coupling(EPC) between the 
TT and 7r*bands is responsible for the peculiar properties observed in the experiments^S'^'^. 

The dimensionless electron-phonon coupling for a mode u at momentum q is given by: 



where Ug^ is the phonon frequency, N„{ef) is the density of states per spin channel at the 
Fermi level, and i and j are band indices. The term 



is the electron-phonon coupling matrix element that describes the scattering of an electron 
from band i to band j due to the phonon u with wavevector q. The quantity Aqj, depends on 
the doping of the system through the shift of the Fermi level and the subsequent change in 
Na-{ef). Furthermore, it (weakly) depends on the doping through the variation of g^^ (k+qij- 
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In the following we are interested in the contribution coming from the matrix elements, 
therefore, in order to lift the dependence on the phonon frequency, we calculate directly 
(k + q, j|AV^,^|k, i). In particular, we will concentrate on the coupling of the tt and vr* 
bands with the highest optical phonon branch (HOB) at F {E2g mode) and at K (A'^ mode). 
We define 

7r,7r* 

{Dl) = J2\{K,J\AVrE,JK,^)\'/A (3) 

and 

(^K)r* = |(2K + k,7r*|AyKA'jK + k,7r)|V2, (4) 

where the sums are performed over the two times degenerate vr bands at K. In the limit of 
zero doping, (-Dp) is equal to {Dy)f as defined in Ref. 
k — > 0, (-D^)]^''* is equal to {D^)f as defined in Ref. 
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27. In the limit of zero doping and 



In fact, for small k the matrix 



elements between tt and tt (or between tt* and vr*) are zero (see note [24] of Ref. 1231). In our 
earlier publications^^ we called these quantities electron-phonon coupling. In view of the 
above definitions (Eqs. ([T]) and ([2])) it is more precise to call (-Dp) {^k)^'* "square 
of the deformation potential" (or, to be very precise, the "average squared deformation 
potential of the vr-bands"). We focus our attention on (-Dp) and (-Dk)^'^*) because these 
two quantities are the ones responsible for the intensity and position of the peaks in Raman 
spectroscopy^, the kinks in ARPES and the phonon slope close to the Kohn-anomalies^^. 

In graphene and carbon nanotubes, the deformation potential has been usually obtained 
from tight-binding Hamiltonians computing the change in the nearest-neighbor hopping en- 
ergy due to a lattice distortion^'^. In this approach the deformation potential, has been 
always considered a constant with respect to the electron or hole density. This approxima- 
tion, although not justified microscopically, is widely \ised!^>^^^^^^^. Moreover even 
ab-initio calculations using Density Functional Theory (DFT) in the local density approx- 
imation (LDA) apparently confirmed that the deformation potential is weakly dependent 
on doping. However, including effects of electron-electron correlation, we will show that 
for the K (A'^ mode) mode, this approximation breaks down and (-DK)k'^* "^^^ change by 
more than 40% just varying the electronic distribution (i.e. gated single and multi-layer 
graphene). This fact can be directly probed measuring the Raman D-peak dispersion of 
graphene versus doping as we will discuss in the following. 

Recently it has been proven^>^ that in (neutral) graphene {D^)^^q is strongly affected 
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by electron correlation. In our previous work—, we have shown that DFT-LDA or DFT- 
GGA underestimates {DI^)^^q by almost a factor of two. The electron-electron correlation 
can be included on the level of the GW-approximation obtaining a deformation potential 
which reproduces the Raman D-line dispersion and the phonon slope around K within few 
percent*. On the contrary, the deformation potential of the T-E2g mode was shown to depend 
very little on electron-electron correlation^''. In this letter, we use the GW-approximation 
to calculate the variation of the deformation potential with doping. 

Doped graphene can be created in single layer field-effect transistor (FET) based ex- 
periments where an electron concentration up to 3 ■ 10^^ cm~^ electron can be realized, 
while higher dopings are obtained using intercaled graphite^^. In order to simulate doped 
graphene we employed a slab-geometry, i.e., bulk geometry with large distance between the 
layers, changing the number of electrons in the unit cell and then compensating the neg- 
ative/positive charge with a uniform positive/negative background, see Ref. [SOj for details. 



The electronic and phonon structure of graphene at different doping levels were computed 
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using DFT-LDA. The deformation potential is obtained using the scheme proposed in Ref. 
based on a frozen-phonon approach by looking at the modification of the electronic structure 
upon displacement of the atoms following a given normal mode. The major advantage of 
this approach is that it can be used with electronic structure methods other than DFT. 
When DFT is used, this approach gives the same result of density functional perturbation 
theory(DFPT). First of all, we investigated the effect of the change in the lattice constant 
a, induced by the doping, on the deformation potential. Using the functional dependence of 
a versus the electron concentration from Eq. (2) of Ref. |3ll, we calculated the deformation 
potential for different doping level with the corresponding lattice parameters. In panel b of 
Fig. [1] we compare the deformation potential calculated with and without lattice relaxation 
for different electron/hole doping. The difference between the two results is small when 
compared with renormalization effects that we are going to describe below. Therefore, in 
order to make the analysis simpler, we performed all the calculations with the graphene 
experimental lattice constant. 

On the level of the LDA, the variation of (D^^)'^^' and (-Dp) with doping is very small. 
In the following, we will introduce correlation effects beyond DFT-LDA on the deformation 
potential. We included these effects in GW approximation'^'^ that has been successfully 
applied in the study of graphene and graphite^"^ and its compounds^"^. First of all we 
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FIG. 1: Panel o: squared deformation potential for the K-A'^ phonon between the vr bands {D'^)^^q 
in different approximations. Circles are GW results; squares are GW{0) ones, keeping the screened 
potential Ty(0) fixed to the undoped case; triangles are B3LYP results. Panel b: variation of 
(-^K)k=o LDA approximation, including or not lattice relaxation (triangles and pentagons, 
respectively) . 
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studied the quasiparticle band-structure versus doping (see also Refs. 
be subsequently used to calculate the D-peak dispersion. In table [H we report the change in 
the gap between the tt and tt* at the high symmetry point M, AeM, as function of doping. 
This quantity is directly related to the optical properties of graphite and graphene-based 
materials, and provides an alternative way to measure quasi-particle renormalization effects. 
The strongest renormalization effect (compared to the LDA-gap) is present for zero doping. 
Electron/hole doping rapidly decreases the GW renormalization of the quasi-particle band 
structure. The same is true for the Fermi velocity vp- Then from the change in the quasi- 
particle band structure upon atomic displacement, we calculated (-Dk)^'^* same way 



as was done for LDA (see ref. 



21). We found that for {Dl) the GW result of ref 
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mainly unaffected by the doping level^^. The situation is completely different for (-D^) 
In Fig. [H and table [H we report the value of (-DK)k=o a function of the doping in different 
approximations. In LDA (panel b) it is almost a constant, but the situation is completely 
different at the GW level (panel a): (-DK)k=o increased more than 80% with respect to the 
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TABLE I: Squared deformation potential for the K-A'^ phonon versus doping, in GW approxima- 
tion. (£'K)kIo and AeM are defined in the text, the LDA value of AeM at zero doping is 4.01115eF. 
D-slope is the slope of the Raman D-peak dispersion (peak energy versus laser frequency), see the 
text. 
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LDA result (at zero doping) and it acquires a strong doping dependence. The renormalized 
deformation potential rapidly decreases with electron or hole doping and gets close to the 
LDA value at large doping. 

In the following, we discuss the origin of the strong doping dependence of (-DK)k=o within 
the GW-approximation. Both the Green's function, G, and the screened Coulomb poten- 
tial, W, are doping dependent. In order to disentangle the two effects, wc performed test 
calculations within the GW approximation keeping the screened Coulomb interaction fixed 
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to its value at zero doping. The result is shown in Fig. [T] The doping dependence of the 
deformation potential is reduced by about a factor of two. Therefore, we conclude that the 
role of the screening is comparable to the correction coming from the Green's function vari- 
ation. In fact the shift of the Fermi level due to doping, leads to a suppression of transitions 
from 71 to TT* states. This affects the deformation potential both through a change of the 
screening and of the Green's function. 
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FIG. 2: Electronic gap at M point (Aem) as function of doping in different approximations: GW 
with fixed screened Coulomb interaction, GW, and LDA 



We compare the result for deformation potential with the one of the AeM gap- For the 
latter we found that all the doping dependence is due to the screened Coulomb potential 
and not due to the change in the Green's function, see fig. [21 while in the first both effects 
contribute to the its behavior with doping. This is due to the fact that, for the K.-A[ 
phonon, GW not only renormalizes the bare Green function lines but also introduces vertex 
corrections that behave different with the doping. However in general this is not true for 
other phonon modes, for instance for the T-E2g mode vertex corrections have been proven 
to be negligible^^. 

At this point it is instructive to check the performance of DFT with hybrid functionals. 
We have performed calculations with the B3LYP functional^i^. Apart from an overestima- 
tion of (-D|j-)j^^Q by about 25% (which could be corrected for by diminishing the percentage 
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of Hartree-Fock exchange in the functional), the calculation reflects rather the doping de- 
pendence of the GW-calculation with constant screening than the dependence of the full 
GW- calculation. This can be understood, because the B3LYP functional contains screening 
on a very simplified level not suitable to describe extended systems and metals. 
Now we want show how the variation of the deformation potential affects the Raman D 
peak dispersion and the splitting of the 2D line in multilayer graphene. The dispersion 
(peak energy versus laser frequency) of the Raman D and 2D lines in graphene is conve- 
niently described by the double-resonant Raman model^i. In order to calculate the phonon 
dispersion we need the deformation potential for a phonon wavevector K-|-q {{Dl^^^)^^*). 
In our earlier works we assumed this as a constant in q and k. The calculation for a finite q 
requires the use of very large supercells, which are challenging for the GW approach. On the 
contrary, we can test the dependence of the deformation potential on the electron wavevec- 
tor k. Indeed, we calculated (D^)'^^ for k varying along the line T — M". We found that 
(-^K)k'^* is almost constant in the doped case and it is slightly varying in the undoped one 
(within 10%), the result is reported in fig. [3l This finding justifies our previous assumptiort^I 
of a constant square deformation potential that was also verified by direct comparison with 
the experiments*. In the double-resonant Raman model^^ the D-line dispersion is propor- 
tional to the phonon slope around K and thus proportional to (D^^)^^* . Furthermore, it 
is inversely proportional to the slope of the tc/tt* bands and thus inversely proportional 
to Aem- In our previous work using the result at zero doping, we were able to reproduce 
completely ab-initio the Raman D-peak— dispersion (peak position versus laser frequency). 
With the information on the doping dependence of {Dl^)'^^* and of Ae^, we have calculated 
the Raman D-peak dispersion as a function of the doping with the approach described in 



ref. 



271 . In Fig. m we report the resulting slope of the D-peak dispersion obtained as a linear 



fit of the dispersion between 1.0 and 3.2 eV laser energy. The D-peak dispersion is almost 
symmetric with respect to electron/hole doping, and it has its maximum at zero doping. 
Due to this strong variation with doping it can be used also to detect experimentally the 
charge state of a graphene sample. Finally we suggest a simple way to verify our result in 
multi-layer graphene using a single laser—. In multilayer graphene the 2D peak splits in 
different sub-peaks and this splitting is proportional to the D-peak dispersion with the laser 
frequencyiiii^. Therefore we expect that measurements of the 2D peak splitting as function 
of doping can highlight the strong variation of the D-peak dispersion, due to the squared 
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FIG. 3: Square of the deformation potential, {D^)^ , with k varying along the line T — M, 
perpendicular to T — K. Four cases are reported: circles and up triangles are GW and LDA 
results for the undoped case; down-triangles and rhombus are GW and LDA for high doping 
An = 9.5410^^cm^^. The upper labels refere to the corresponding electronic gap (the energy 
difference between vr* and n bands at K+k according to GW calculations at zero doping). 

deformation potential, as predicted in this work. 

In conclusion, we have shown that in graphene {D'^)'^^ can be easily tuned with doping. 
This means that the electron and phonon interaction cannot be described by a simplified 
Hamiltonian with a fixed deformation potential as it is commonly done. Our findings can 
be easily verified experimentally by measuring the doping dependence of the Raman D- 
peak dispersion, of the splitting of the 2D peak in multi-layer graphene and of the slope 
of the highest phonon branch close to K. The present findings open the possibility to use 
the slope of the D-peak dispersion as a simple probe for the electrons/holes doping in 
graphene. The fact that the deformation potential is not a constant has implications also 
in the realization of graphene-based electronic devices where it is known that one of the 
limitation for ballistic transport is just the coupling between electrons and phonons^i"-. 
Tuning electron-phonon coupling by doping could boost electronic properties of graphene 
based devices. Concerning the puzzling discussion on the size of the EPC in graphene, our 
result puts another piece in support that EPC has to be larger than the LDA one and doping 
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FIG. 4: Change of slope of the Raman D-peak dispersion versus doping . Arrows indicate the 
equivalent doping level for the KCg, KC24, KC36 intercalated graphite. 

dependent. Finally, this new way to tune electron-phonon coupling can play a fundamental 
role in all the experiments where physical phenomena are directly related to doping as for 
instance the phonon renormalization^, phonon linewitdh^, or the radiative decay of excitons 
in nanotubes^. 
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